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Abstract 

Given a sequence of integers a,j , j > 1 , a multiset is a combinatorial 
object composed of unordered components, such that there are exactly 
aj one-component multisets of size j. When aj X j'^ 1 y^ for some 
r > 0, y > 1, then the multiset is called expansive. Let c n be the 
number of multisets of total size n. Using a probabilistic approach, we 
prove for expansive multisets that c n /c n +i — > 1 and that c n /c n +i < 1 
for large enough n. This allows us to prove Monadic Second Order 
Limit Laws for expansive multisets. The above results are extended 
to a class of expansive multisets with oscillation. 

Moreover, under the condition aj = Kj r ~ 1 yi + 0{y y i), where K > 
0, r > 0, y > 1, v G (0, 1), we find an explicit asymptotic formula for 
c n . In a similar way we study the asymptotic behavior of selections 
which are defined as combinatorial objects composed of unordered 
components of distinct sizes. 
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1 Summary and Historical remarks 



Given a sequence of integers a,j > 0, j > 1, a multiset is a combinatorial 
object of finite total size composed of unordered indecomposable components 
such that there are exactly a,j single component multisets of size j. There is no 
restriction on the number of times a component may appear in the multiset. In 
view of this, let 

fi n = j?? = (?7i , ?7 2 , . . . , r) n ) : ^2 Mi = n and Vj > for all j J 

be the set of unordered integer partitions of an integer n. Any multiset of 
total size n has a component count vector if = (771,772, • • • ,Vn) contained in 
Q n . Here rjj = rfj^ is the number of components of size j in the multiset of 
size n considered. (For more details about multisets see e.g. [2, 9]). 

Let c n be the number of multisets of size n determined by the above 
parameters a,, j = 1, . . . ,n. We will prove an analytic identity which will 
be used to extract information about the growth of c n , as n — > 00. It follows 
from the definition of a multiset that the number of multisets with a given 
component count vector 77 = (771, 772, ... , i] n ) G f2 n is 




Hence, the number of multisets of size n has the Euler type generating func- 
tion g: 

00 

n=0 

00 Tl / x 

n=lrjen n j=l \ / 

00 

= 11(1-^. N<!" (!) 

We observe that combinatorial objects that are given by the generating func- 
tion (1) are also called weighted partitions ([8]). 
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The truncated generating function 

n 

9n(x)=ll(l-^y ai , \x\<l (2) 

3=1 

has Taylor expansion g n (x) = YlT=o c k,nX k , where Ck, n = c n for k < n. For 
a fixed n, the series expansion of the function g n (x) converges for all x G C, 
\x\ < I. We set 

x = x ( a , a) = e - a+2ma 
for some real numbers a and a. Then we have 

jT g n (x)e- 2man da = £ ^ c^e^ 2 ™^-^ da = c n e~ n % (3) 



where we have used 



i 

e 2niam da 



o 



1, if m = 0, 

0, if m ^ 0, m G Z. 



Substituting (2) into (3) gives the desired identity w.r.t. the free parameter 
a ER: 

/l n 
II ( X ~ e~ ja e 2nia:i ) ~ aj e- 2nian da. (4) 

A more probabilistic interpretation of (4) can be given. We have 

1 n / i p—ja \ a i 

! ' > -2nian 



e-'" m da 



n -i 

= e na \\{l-e~ ja y a3 / 0(a) e" 2 ™ n rfa, (5) 
.7=1 Jo 



where 0(a) is given by 



0(«) =n^ 

^(a), a G R, (6) 
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for functions <pj defined by 

/ 1 _ e -3<r \ a J 

M°)={ 1 _ e _ jl , eMaj ) , (7) 

Using the combinatorial identity in Lemma 2.46 in [9] it is easy to see that 
for any a > 0, 0j is the characteristic function of a random variable Xj given 
by 

HXj = Jl) = (° J + / ~ (1 - e-'T e~ lj °, 1 = 0,1,2,.... (8) 

Consequently, 0(a) is the characteristic function of Y n := ^2™ =1 Xj, where 
the Xj, j = 1, . . . , n are assumed to be independent. Therefore, 



i 

2nian 



<j)(a)e- 2nian da = F(Y n = n). (9) 



Combining (9) with (5) we arrive at the desired representation of c n , 
which is in the core of the probabilistic method suggested by Khintchine in 
1950-s ( [20], Chs IV, V) for asymptotic enumeration in the urn models of 
statistical mechanics. The history related to the method is discussed in [15]. 
We note that Khintchine- type representations were subsequently rediscovered 
in independent ways by many authors in a variety of seemingly unrelated 
contexts. In particular, observe that (5) can be derived from equation (134) of 
[3], the latter being based on the conditioning relation (see [2]). In conclusion, 
the recent paper [22] should be mentioned which treats probabilistic methods 
for enumeration as transforms of generic random variables (in our setting Zj) 
into specially constructed independent random variables (in our setting Xj). 

It follows from (8) that the r.v. j~ 1 Xj is negative binomially distributed 
with parameters a, and e _<7J , a > 0. This produces the formula for the 
expectation of the sum Y n : 



n . _ ja 



M n :=EY n = J2^-J,- (10) 
i=i 

Further on, except for Lemma 1 in Section 2, we will assume that a — a n > 
is chosen in the unique way so that 

M n = n. (11) 
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The fact that a can be chosen in such a way follows from observing that M n 
decreases monotonically from oo to as a ranges from to oo, while n is 
fixed. The variance of Y n is 

B; : =Va r( y n ) = E n_ a e-^V <12) 

j=\ 

We will study the probability in (9) analytically and thereby obtain in- 
formation about asymptotic behavior of c n , as n — > oo. 

It is natural to suppose that, under some asymptotic conditions on the 
parameters aj, j > 1, a local limit law should hold for Y n in (9). Asymptotic 
enumeration of multisets using this approach was apparently initiated in [13], 
which was preceded by works of Freiman ( see [24]) on the development of 
Khintchine's method. Note that asymptotics of different statistics related to 
integer partitions (the case aj = 1, j > 1) was studied by many authors(see 
e.g.[16, 23]). 

In this paper we will initially assume that 

aj x J r ~y , j ^ oo, r > 0, y>l, (13) 

where we write aj x bj for sequences aj, bj when there exist constants 
D 1 ,D 2 > such that D 1 < aj/bj < D 2 for all j > 1. Although for y > 1, 
the rate of growth of aj is exponential (but not polynomial) such multisets 
will be called, following [5] , expansive. This is in view of Bell-Burris Lemma 
5.2 in [7] which tells us that for y > 1, the asymptotic behavior of the num- 
ber of multisets with aj ~ j r ~ 1 y-' is the same as that of assemblies with 
aj ~ f~ x ■ (Here and in what follows a n ~ b n for sequences a n , b n means 
that lim ?woo a n /& n = 1). 

Provided the parameters aj satisfy (13), we will prove that the normal 
limiting law for the sum Y n holds, leading to f(Y n = n) ~ (27r£^) -1 / 2 in (9). 
This local limit law, the definition of a n , and (5), will imply our Theorem 1 
below. A result implying Theorem 1 when in (13) y = 1 and x is replaced 
by ~, was obtained by Richmond in [25] and applied for enumeration of 
partitions of n into primes or powers of primes. The first probabilistic proof 
of the Hardy- Ramanuj an formula for partitions ( the case aj — 1, j > 1) was 
given by Freiman in 1950-s (see Section 2.7 of [24].) Quite recently, a different 
probabilistic proof of this formula was published in [4]. A comprehensive 
study of the asymptotics of integer partitions was made in [?]. 
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Theorem 1 appears to be new for y > 1. Note that throughout the paper 
we assume, unless it is said otherwise, that all asymptotic expressions are 
valid under n — > oo. 

Theorem 1 Assume that (13) holds. Then the number of multisets is asymp- 
totically 

where a n is given by (11). 

We now formulate an extension of Theorem 1 to a wider than (13) class 
of parameters aj. Corollary 1 below is an analog of Corollary 1 in [14] for 
expansive assemblies. 

We write qi(n) z< •(n) ^ ^(n), if there exist positive constants 71, 72, 
such that 7i<7i(n) < m(n) < 72<?2( n ), n > 1 • For given < r\ < r 2 and y > 1 
define the set jF(ri, r 2 , n) of parameter functions a = dj, j > 1 satisfying the 
condition 

rv^i^rv, 2/ > 1, j > 1. 

Corollary 1 For an arbitrary r > 0, < e < r/3 and y > 1, the conclusion 
of Theorem 1 is valid for all parameter functions a G jF(2r/3 + e, r, y). 

It is interesting to observe that in the case y — 1, our condition a G 
jF(2r/3 + e,r,y) implies the condition (i), p. 1084 of Richmond [25]. This 
latter condition is sufficient for the asymptotic formula for partitions of n into 
primes developed in [25]. Generally speaking, multisets with a G J-(ri,r2,y) 
may be called "expansive with oscillation" . 

Theorem 1 and Corollary 1 are proved in Section 2. In Section 2 we also 
derive asymptotic estimates for a n and that are used in (14). 

A multiset satisfies a monadic second order logical limit law if the proba- 
bility that a random representative of the multiset of size n satisfies a monadic 
second order sentence converges, as n — > 00. Compton [10, 11] showed that to 
prove that a class of relational structures such as multisets satisfies a monadic 
second order limit law, it suffices to know about the growth of the number 
of structures c n of size n. The next corollary from Compton's theorem was 
used in [7] to prove logical limit laws. 
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Theorem 2 [Corollary 8.1 of [7]] Suppose that aj are the parameters of 
a multiset such that 



~ y for some y > 1. 

c n+l 

If y > 1, then suppose further that there exists N > such that 

—— < y 1 whenever n > N. 

Cn+l 

Then the multiset has a monadic second order logical limit law. 

Based on a Tauberian theorem of Schur, Bell and Burris ([7], Theorems 9.1 
and 9.3) and Bell [6] derived general sufficient conditions on the asymptotics 
of aj which imply the hypotheses of Theorem 2. Note that the condition on 
the a.,- obtained in [7] covers the particular case aj ~ j r ~ 1 y^, y > 1 of (13). 
Combining Theorem 1 with the asymptotic estimates in Section 2, we prove 
in Theorem 3 below the validity of the conditions of Theorem 2, in the case 
(13) that is not covered by the Bell-Burris sufficient conditions. 

Theorem 3 Suppose that aj is a sequence obeying the condition (13). Then 
the corresponding multiset has a monadic second order logical limit law. 

A result similar to Theorem 3 is obtained for logarithmic structures in 

[?]■ 

Moreover, we are able to weaken the condition (13) of Theorem 3: 

Corollary 2 For all parameter functions a e JF(2r/3 + e,r,y), where r > 
0, y > 1 and < e < r/3, 

= y- 1 exp (-8 n + o (5 n )) , 

Cn+l 

where 5 n = a n — logy — > 0. 

In particular, Corollary 2 implies that c n /c n+ i — > y > 1, for a e 
jF(2r/3 + e,r, 1). A similar result was shown in [14, 15] for certain reversible 
coagulation-fragmentation processes. From an analytical point of view the 
latter processes are equivalent to expansive assemblies (see [15]). 
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In view of the above definition, we may consider multisets as unlabelled 
decomposable combinatorial structures. We call labelled decomposable com- 
binatorial objects assemblies, a term used in [3], see also [2]. Sufficient con- 
ditions to have monadic second order logical limit laws were given for both 
multisets and assemblies in Theorem 6.6 of [10]. Let rrij be the number of 
labelled components of size j and let aj = ^f. For assemblies the basic 
representation (5) becomes for an arbitrary choice of a 



see (2.24) of [14] or (125) of [3]. The method of proof of Theorem 3 and the 
comment in the last paragraph gives monadic second order logical limit laws 
for assemblies whenever a G jF(2r/3 + e, r, 1), r > 0. 

Theorem 3 and Corollary 2 are proved in Section 3. 

The problem of "factorisatio numerorum" can be put in the framework of 
enumeration of multisets. The following description of factorisatio numero- 
rum is taken from [21]. An ( additive) arithmetical semigroup is a free commu- 
tative semigroup G with identity element 1, generated by a countable set P 
of "prime" elements, and equipped with an integer- valued "degree" mapping 
d such that 

(i) <9(1) = 0, d(p) > for all p G P. 

(ii) d{ab) = d{a) + 8(b) for all a,b G G. 

(Hi) The number G # (j) of primes of degree j in G is finite for all integers 



Multisets can be put into the framework of arithmetical semigroups by 
letting the operator d stand for the size of the multiset and defining the 
product of two multisets to be their disjoint union. The identity element 1 
is then just the empty multiset with total size 0. 




where 




r > 0; 



J- 
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Let f(n) be the total number of unordered factorizations of elements 
g G G with d(g) = n. Then [21] shows that 

oo oo 

J2f{n)x n = Y[{l-x j )- G * ij \ \x\<l. 

n=0 j=l 

This equation is just (1), except that c n has been replaced by f(n) and cij by 
G*(j). A typical example considered in [21] is polynomials over finite fields, 
for which = q* for some prime q > 1. 

We are able to extend the results of [21] and give asymptotic results 
for "factorisatio numerorum" when aj = = Kj r ~ 1 yi + 0(y v i) for 

K > 0, v G (0, 1) and r > 0, y > 1. This involves getting precise enough 
estimates of a n in order to derive first order asymptotics of c n . We restrict to 
the case y > 1 , as then a fairly simple argument using the Poisson summation 
formula is effective. 

Theorem 4 Assume that aj = Kj r ~ 1 y j + 0(y uj ), where K > 0,r > ; 
y > 1, and v G (0, 1). Then c n has asymptotics 

c n ~ /€iy n n-( r+2 )/ 2 ( r+1 ) exp {K 2 n r '^) 

for positive constants K\, k 2 . Moreover, 

K2 = r -±l(KT(r + l)) l K r+1 \ (15) 
r 

where Y is the gamma function. 

For r = 1, (15) recovers the asymptotic formula in [21]. 
Theorem 4 is proved in Section 4. 

Remark 1 It is known ([9], p. 34) that the generating function g for c n can 
be written as 

g(x) = exp S** (re), \x\ < 1, (16) 

where 

oo 

S*{x) = Y,a*x j , \x\<l 
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is the so-called star transformation of the generating function 



S( x ) = ^^ttj^i \ x \ — i 
3=1 

for dj, namely 

ik=j 

The representation (16) says (see e.g. [15]) that g can be viewed also as 

m* 

a generating function for the parameters a*j — -jf , j > 1 of assemblies. By 
Lemma 5.2 in [7], we have that the asymptotic formula in Theorem 4 for 
enumeration of expansive multisets with y > 1 is also valid for enumeration 
of assemblies with the same parameters aj. In this connection observe that, 
under the assumption (13), the orders of the quantities 5 n , B\, pi(n) found in 
Section 2 appear to be the same as the ones in [15] for expansive assemblies. 
Summing this up, we see that the asymptotic behavior of expansive assemblies 
and multisets is alike. We will show further on that the same is true also for 
selections. 

We define the selections determined by the parameters aj to be those 
multisets for which no component type appears more than once. For example, 
if aj = 1 for all j then a selection is an integer partition with distinct parts. 
Let c n denote the number of selections of size n determined by the aj. Then 
the generating function g for the c n is 

oo 

~g(x)=l[(l + xlp, \x\<l 

3=1 

and analysis similar to the one for multisets gives that in this case j~ l Xj is 
a binomial r.v. with parameters a,j and 1 ^p"(-jCT) ' w ^ ere °" > is arbitrary. 
Consequently, we have 

Theorem 5 Assume that aj satisfy (13). Let a n be chosen in such a way 
that 

n . 
J a i e 

M n := > ^ — = n. 

3=1 
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and define B n by 



J a j e 

3=1 y ' 
Then the number of selections is asymptotically 



£n~ i — \l{l+e-^) a \ (18) 



Moreover, if we assume that aj is as in Theorem 4 then c n has the same 
asymptotics as c n , with a different constant k\. 

We sketch the proof of Theorem 5, which is similar to the proof of The- 
orems 1 and Theorem 4, in Section 5. 

The classic example of an expansive multiset is integer partitions. For 
partitions of an integer, aj = 1 for all j, so that r — y — 1. We can derive 
the Hardy- Ramanuj an formula giving asymptotics of c n for partitions from 
Theorem 1 by using well known properties of the Euler generating function 

F ( x ) : = UT=i (! - &T 1 , M < L Since Jo°° f^ du = T ( see t 18 ]' Formula 
3.411-7), we apply the Euler-Maclaurin formula (described in detail in [17]) 
to obtain from (11) 

~- i 

Consequently, 



TT 2 


1 


6^ 


2(7,, 


7T 


1 




4n 



)(n" 3 / 2 ). (19) 
The equality (8.6.1) in Section 8.6 of [19] gives 

F(e -~) (20) 
Finally, Theorem 1 and (20) produce 

°n ~ /^— FTT U~ exp — . (21) 



V ^B 2 V27r7 ' ^ V6^n 
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An asymptotic analysis using (19) which we do not present here shows that 
the asymptotic relation (21) still holds when a n is replaced by 7c/^/6n. Fur- 
thermore, 



3 D 2 r u2e ~ u a -k 2 

n n J (l-e-«) 2 3 



(see [18], Formula 3.423-3). The above analyis results in the Hardy- Ramanuj an 
formula 



e Cy/n 



with C = tt(2/3) 1/2 . 



Examples of expansive multisets, most of them with r = 1, can be found 
in [7], [9]. The simplest example is the class of finite /c-colored linear forests, 
which has aj = & , so that r = 1, y = k. We give an example with r = 1/2, 
y = 2. 

Example 1 Consider the linear forests in which every tree on j vertices is 
2-colored with colors red and blue in such a way that it has exactly [j/2] red 
vertices and j — [j/2] blue vertices. Then aj = (yj 2 \) ~ \/2/7r2 : 'j~ 1 / 2 . 

We may generalize the last example to get multisets with any r£ [1/2, 1] 
as follows. 

Example 2 Consider the forests in which every component of size j is com- 
posed of a linear tree on [j a \ , a e [0, 1] vertices and a cycle on the j — [j a \ + 1 
vertices, in such a way that one end vertex of the tree is identified with one 
vertex of the cycle. Call such components (= graphs) lollipops. Suppose we 
are considering 2-colored lollipops, such that the number of red vertices in 
the tree is [j a /2] and the number of blue vertices is j a — [j a /2]. There is no 
restriction on the number of blue/red vertices in the cycle. The number of 
2-colored lollipops is 

This example has r = 1 — a/2, y — 2. 

The next example is a natural case where the multiset satisfies Theorem 3 
but not the conditions in [6, 7]. 
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Example 3 Consider linear forests which are k colored. If a tree is j vertices 
long and j is even, then it may be k-colored in all & possible ways. If j is 
odd, then the first vertex is always red and the remaining j — 1 vertices may 
be colored in all W~ l possible ways. Then a,j = W if j is even and aj = W~ l 
if j is odd. Therefore, (13) holds with r — l,y — k, D ± — k~ 1 ,D 2 = 1. 

Example 4 Finally, note that r = 2, y = 1 corresponds to plane partitions; 
see [1]. 



2 Asymptotics for expansive multisets 

In this section we prove Theorem 1 and Corollary 1. Recall that we assume 
here, with the exception of Corollary 1, that a,j obey the condition (13). We 
first derive an expansion for the characteristic function given by (6) to 
general precision. For any a > 0, we define the quantities pi = pi(n) for 
I > 3 by 

n oo 

P^EA-E^ 1 ^- (22) 
j=l k=l 

Lemma 1 For a fixed n and any integer s > 3, the function <fi can be ex- 
panded as 



<f)(a) = exp y2mM n a - 2n 2 B 2 n a 2 + ^ ( 2 ^) Pi a i + ( a s p s ) J , a -> 0, 

where M n and B n are given by (10) and (12). 
Proof The definition (7) implies that for all a G R 

™ / 1 _ -ja \ a 3 
3=1 V 7 

= exp (j^ % ( ^g (1 " e-^) - log (1 - e-^") ) j . 



13 



The logarithms may be expanded in Taylor series as a > and a G R are 
fixed, giving 



0(a) = exp + ^ 

0=1 V k=l k=l 



We make use of the Taylor expansion with s > 3 
e 2 "*" = 1 + 2nijka - 2 I r 2 /fcV + £ (2 "^ a) + O(fk'a-), a^O, 



which holds uniformly for all j, k > 1, to get 



/ e -jka 



0(a) = exp 5Z a M S 



k 

3=1 \k=i 



2nijka - 2n 2 3 2 k 2 a 2 + ]T + ( a »^ 



Z=3 



/ . . 7a 9 e ja n 9 7 2 a,e J,cr 9 £-4 (2m) 1 pi , _ 
= exp 2tu> -J-l — — a -2-K 2 y , 3 ■ a 1 + 0(p s a s 

j=l 3=1 V 7 1=3 



= exp ( 2™M n a - 2n 2 B 2 n a 2 + ]T ^^a z 

V 1=3 L 



+ 0(a a p 9 ) , a ^ 0. 



In what follows we set a = a n determined by (11) and define 5 n by 
$n '■= °~n — logy, V > 1- In proving Theorem 1 we will apply Lemma 1 and 
so we need estimates of 5 n , B n and p 3 . (We use the ability to expand 0(a) to 
higher order precision than s = 3 for the proof of Theorem 3 in Section 3.) 

Lemma 2 5 n x n" 1 /^ 1 ), B 2 n x n ^^ r+1 \ and Pl (n) x n ^^ r+1 \ r > 0, 
for all I > 3. Moreover, there exists N > suc/i i/iai 5 n+ i < 5 n; -B 2 < B 2 +1 , 
and pi{n) < pi(n + 1) whenever n > N. 

Proof We first prove some preliminary facts about 5 n . Let Di,D 2 > be 
constants such that D\f~ x yi < a 3 - < D 2 j r ~ 1 y 3 \ j > 1, y > l,quadr > 
0. Since <r n > 0, n > 1, by (11), we deduce from (10), (11) that n > 
Y^j=i3 a 3 e ~ 3 ' Tn > Di YTj=i j r e~ j5n , implying that 5 n > for n large enough. 
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Suppose that there exists a constant e > and a subsequence n k — > oo such 
that 5 nk > e. Then, again by (10), (11), 

D nk 

n k < - \ — - y^fe' 36 = 0(1), n k -> oo. 

1 — v~ 1 e~ e ^— ' 
y i=i 

We therefore must have 5 n — > 0. Next, we derive from the inequality 

n n 

n > D l J2f e ~ j5n > Die~ n5n ^ J r , r > 

that n5 n — > oo as n — > oo. 

With the help of these facts we further get 

n<D 2 y^ <5 n r 1 D 2 J2 K ' _ lS S n ~S n r X D 2 / -dx. 

3=1 y 3=1 ° 

Since the last integral is bounded, we conclude that S n < D A n~ l ^ r+l \n > 
1, where D 4 > is a constant. On the other hand, 

n> D 1 J2f e ~ j5n ~ / a: r e"^x. 

This gives <5 n > D 3 n _1 /( r+1 ), n > 1, where -D 3 > is a constant. We have 
shown that 5 n x 

For any I > s > 0, arguments similar to those above show that 

Y / a l e . x y 3 x 5 -r-l x n (r+l)/(r + l) (24) 

' (1 — e -3^njS Z-^ (I _ y-3 e ~3°ny n v 7 

The last asymptotic applied to (12) results in the stated asymptotics for B\. 

We now show that pi, I > 3 has the same asymptotics as in (24). We 
have 

n oo 

Pl x y^jr+i-iy^ k i-i e -5 n jk y j-jk 

3=1 k=l 

~ 5 n r "' X] (^" r ~ 1 / ^ +, - 1 !/^ (1 - fc) e-*^) , r > 0. 
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The integral in the last expression is < T(r + I) for all k > 1, with equality 
for k — 1, since n5 n — > oo, <5„ — > 0, n — > oo, while i/*' 1- '' < 1, k > 1, y > 
l,z > 0. Thus, the last series in (24) converges which implies 

^ x x (25) 

Lastly we will prove that <r n , and so 5 n , is eventually monotone decreasing 
in n. Suppose that cr n +i > °n for some n. Then, it would follow that 

< 



X — g~jo"n+l \ — g— jCTn 

for all j > 1, and consequently 

^ ja je -^ < jaje-^" = ^ (n + lK+iy^e^^ 
i=i j=i y 

In view of the established asymptotics for 5 n , the last term in the preced- 
ing inequality tends to for all y > 1. We therefore must have <7 n+ i < a n for 
sufficiently large n which implies the same inequality for 5 n . 

The derivative ^ [e~ x /(l — e~ x ) 2 } is negative for x > 0, so that if n is 
sufficiently large, 



2 

n+1 



and similarly, pz(n) < pi(n + 1) when n is large. 

Lemma 3 below proves a local limit theorem for the probability in (9). 
Define the sequence a (n) by 

a (n) := ^ r+2)/2 logn x n -(r+2)/2(r+i) logn _ (2 6) 
We will express the integral in (9) as 

r l f l/2 

T = T(n):= / <P(a) e~ 2nian da = / 0(a) e" 2man rfa = T x (n) + T 2 (n), 

JO 7-1/2 

(27) 
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where the middle equality follows from the periodicity of <p(a), as defined by 
(6), (7), and where 



/a (n) 
<j)(a) e~ 2nian da, (28) 
-ao(n) 

/-ao(n) pl/2 
<P(a)e- 2nian da+ / 0(a) e~ 2nian da. (29) 
-1/2 Ja (n) 

Lemma 3 and the representation (5) prove Theorem 1. 
Lemma 3 

7\ ~ (2vr^)- 1/2 (30) 

and for sufficiently large n, 

T 2 <exp(-Clog 2 n), (31) 
/or a constant C > 0, from which it follows that 

v -1/2 



T = / 0(a)e- 2man cfa ~ (2tt5 2 ) 



Proof The proof of this lemma is similar to the proof of Lemma 6 in [14]. 
Using the expansion of Lemma 1 with s = 3 in the definition (28) of T\ and 
observing that, by virtue of Lemma 2, lim^oo a 3 p 3 = and a 2 B 2 — > oo for 
all a G [— cxo(n), a (n)], gives 

rao(n) 

Ti ~ / exp (-2n 2 a 2 Bl) da ~ (27T.B 2 ) 7 . 

./ — ao(n) 

The bound for T 2 starts with the identity for all a G R, 

fr i-«- J 7 .^f-£ai«(i + ^»*y>) 

11 i _ e -j<Tne*maj ^ \ Z-j 2 \ (1- e -^™) / 



Wl = 



All of the logarithms are positive, and log(l + x) > x/(l + c) whenever 
< rr < c for a constant c > 0, so we have 

. ,. .. / a,- / 4e~ jCTn sin 2 (7ra?')\ \ 

Wa)| < exp - JJ i'°g( 1+ (l-e-V ) 
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^2 C 1 a j e~ jan sm 2 (7iaj)\ 

( 4( j n )-i<j<n / 

C 2 f~ 1 e- jSn sm 2 (vraj) , a G R, 

(4a n )- 1 <j<n / 

for some constants Ci,C 2 > 0. In view of the inequality 5 n < 5 n + logy = 
(T n , y > 1, we also have 

|0(a)| < exp j - - C 2 f~ 1 e- j5n sin 2 (vraj) J := exp(-V n (a)) , a G R. 

\ (45„)- 1 <i<n / 

(32) 

Since our 5 n is of the same order as a n in [14] , the argument of Lemma 7 
in [14], gives the desired estimate of 14(a): 

V n (a) x 5~ r x n T / (T+1 ^ > log 2 n, a G [a , 1/2] (33) 

Proof [of Corollary 1] For a G jF(2r/3 + e, r, y), r > 0, < e < r/3, y > 1, 
the arguments in Lemma 2 show that 

n -V(ri+l) ^ ^ ^ n -V(r a +l) > ^-(n+2) ^ £2 ^ ( J-(r 2 +2) > (34) 

£~ (ri+0 =< Pi{n) ■< 5~ (r2+l) for Z > 3. (35) 
We write a = (B n )~ l logn to obtain, as n — > oo, 

agpaH < 73 (log 3 n)e i+2)/2 ^ (r2+3) = 7s(log 3 n)e i - 2r2)/2 - 0, 

since 3r x — 2r 2 > 0. We have shown that the asymptotic (30) is still valid. 
The upper bounds on T 2 are like those in the proof of Theorem 1, with the 
replacement of (33) by 

K ri =< V n (a) ± 5- r > for a G [a , 1/2]. 



< exp — 



< exp 
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3 Logical limit laws for expansive multisets 

Lemma 4 below and the asymptotic 8 n+ i x n^ 1 ^ r+1 ' ) from Lemma 2 show 
that the c n satisfy the hypotheses of Theorem 2 and thereby prove Theorem 
3. 

Lemma 4 If aj x j r ~ 1 y j , with r > and y > 1, then 

° n — y-i e -Sn+o(8 n ) _ 

Cn+l 

Proof We use (5) to get 



Tn-(n+l) CTn+ l TJ / ^ \ ^ _ e -(7l+l)<7„+iy 



Cn+i f = \ V l-e-f" J K ' T{n+iy 

(36) 

where T{n) is defined by (27). Since 

e -(n+l)a n+lan+i x ^ + 1 ^r-l e -(n+l)«„ +lj 

by definition of 8 n , and 

(n+l)5 n+1 x(n + ir^ +1 ), 

by Lemma 2, it follows that 

(i.e-W^i)^^^). (37) 

The second factor in the RHS of (36) may be rewritten as 

We assume that n > N for the iV in Lemma 2, so that in particular <r n +i < <r n . 
Since log(l — x) < —x when x G [0, 1], we have 

" / l-e-^ + i \ a ' < ^ / " e-^-e-J" 
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c e Ja "(jcr n - ja n +i) 
< expl-^a, — 



g-JCTn 

3=1 



= exp I -(a n - a n+1 )^2~ 



3=1 

= exp(-((j n - <j n+ i)M n ) 

where the second inequality results from the fact that e z — 1 > z, z > 
Since log(l — x) > —xj (1 — x) for x G [0, 1], we lower bound (38) by 



J"=l V ' V 3=1 



> exp - > a,- : 

3=1 



= exp I -(a n - a n+l ) ^ 



n+l . , . 



g-J°"n+l 

i=i 

> exp(-(<7 n -<7 n+ i)M n+ i) 

= g^( cr ™ _cr "+i)( n + 1 ) 

Thus, the product of the first two factors of (36) is bounded above and below 
by 

JL / 1 _ g-i^n + i \ a i 

g-0" n < gno- n -(n+l)(T n + l 1 I I \ ^ g~0"n + l _ g-0" n + (o" n -CT n + l) 

- 11 I 1 _ g-ja n J ~ 

3=1 

We bound <7 n — <r n+ i by observing that 
1 = (n+l)-n 



3=1 3=1 

J'=l 
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j 2 cije jrTn 



> in,. n ; .. ; )^ T 
i=i 

for n sufficiently large. 

Thus, recalling that a n — <7 n +i > 0, it follows from (24), applied with 
I = 2 and s — 1, and Lemma 2 that 

a n - a n+1 < O (C 2 ) • (39) 

We have shown that 

-A- / 1 - e~ j<Tn+1 \ flj 

II( !_ e ^ j = e W (-a n + o(5 n )) 

= y' 1 exp (S n + o (5 n )) . (40) 



raa- n -(n+l)cr n+ i 



Because of (36), (37) and (40), the proof will be completed if we show 
that = e°( 5 "\ The definitions (27), (28), (29) along with Lemma 2 and 

Lemma 3 imply that 

T 2(n) 

which gives 

T(n) 1 , 7\(n) -T 1 (n + 1) 



o(5„). (41) 



T(n + 1) 7i(n + l) 

The definition (28) together with (26) and (39) produce 

/a (n) 
|0 n (a)e^ 2 ™ 1 - n+1 («K 2 ™ (n+1) | da + o(5 n ). 
■ao(n) 

Next, Lemma 1, (11) and the monotonicity of pi(n) imply that, for a 
sufficiently large fixed n and a — > 0, 

= exp (-2tt V5 2 + Q s (a, n) + 0{a s p s {n))) 

- exp (-2?r V5 2 +1 + Q s (a, n + 1) + 0(a%(™ + 1))) 
= n (aK 2 ™" 

x 1 1 - exp (-27r 2 a 2 ( J B 2 +1 - 5 2 ) + Q s (a, n + 1) - Q s (a, n) + 0(a s p s (™ + - 
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where we denoted 



1=3 

We now apply (39), Lemma 2 and (24) with / = 3, s = 2 to get 

n < r>2 _ R 2 \ ■/ "j e \ -V "j e 

- n+1 ™ "~" (1 _ e -j<r n+1 \2 (1 _ e -ja n )2 

i=i v 7 j=i v 

^ fa je -^ _. (n + l) 2 a n+1 e^ n+1 ^ 

< 0(C)- (43) 
In a similar way we also have from (22) 

n 00 

0< Pl (n + l)- Pl (n) < {a n -a n+l )Y,J l+1 ^Y, kle ~ 3kan+1 

j=i k=i 

+ {n+ l)2J+r-l e -i„(n+l) 
n 00 

< 0(<Sn)^j ,+r £^~ J ' fc *" +1 

i=i k=i 

n roc 

= 0{8- l )Y,r 1 / ti'e-du 

Vie/?! ieD 2 J « E / 

< O (5~ {l+r) n er ) , />3, e>0. (44) 

Here D x = [1,5"^^] and D 2 = [l,n] \ Z?i. For a G [— a (n), a (n)], (26), 
(43) and Lemma 2 imply 

\(Bl +1 -Bl)a 2 \<O(5 r n +1 \og 2 n)^0. (45) 

Similarly, for a e [— «o( n )? Q; o( n )]; ^ > 3 and all r > 0, 

|(p,(n + l)-p,H)a'| ^(^"VHog^^O, (46) 
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for sufficiently small e > 0. Because of (26) and Lemma 2, it follows that for 
a G [-a (n),a (n)}, 



\a l \pi(n) = O (5 r n (l - 2)/2 log' n) -> 0, I > 3. (47) 

The above discussion reveals the following remarkable feature of the choice 
a G [— a (n), a (n)] in the expansions (42) and in Lemma 1: under this 
choice the terms with s > 3 can be ignored, asn-> oo. Therefore, based on 
the preceding bounds, we get 

\Ma)e- 2man - (j) n+1 (a)e- 2ma ^\ 
< 0{n~ 1 log 2 n) exp (-2vr V5 2 + O (5 r r / 2 log 3 n) ) , 

uniformly for a G [—ao(n),ao(n)]. Now it follows that 
\TAri) -7\(n + l)| ^ 2 , 

and therefore (41) gives T "^"^ = 1 + 0(n _1 log 2 n) + o (5 n ) = 1 + o (S n ) = 
exp (o (5 n ))) proving the lemma. 

Proof [of Corollary 2] We will make use of (34) and (35). The argument 
in the proof of Lemma 2 shows that in the case considered we still have that 
8 n — > as n — > oo and that 5 n decreases monotonically for large enough n. 
Let, as before, r\ = 2r/3 + e, < e < r/3 and r 2 = r. Then, observing 
that the bound in (39) is valid with r replaced by r±, we get that the bound 
in (43) becomes 0(^ 1_r_1 ). Consequently, setting a as before gives (45) 
with r replaced by r\. The left hand sides of (44) and (47) may be bounded 
similarly. 



4 Explicit asymptotic formulae for enumera- 
tion of expansive multisets 

In this section we will prove Theorem 4, which gives first order asymptotics 
for c n when y > 1, K, r > and 

aj = Kf~ 1 y j + 0(y uj ), where v G (0, 1). (48) 
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To approximate a n = logy + S n in the case considered, it is necessary to 
analyze the equation 



7a,-e Ja " 
n - * J 



3=1 



= J2 ja i y 3 e jSn +° J2 ja i y 2je 

3=1 \j=l 
n / n \ 

= Kfe- jSn + o J2i r y~ {1 ~ u)je ~ j5n + 

j=l \J"=1 / 

n 

= ^iffe^'* 5 " +0(1). 
J'=l 

The Poisson summation formula as used in the proof of Lemma 4 of [12] 
shows that for / > — 1, 

n 

J2j l e-^ = T(l + I)*" 1 " 1 +Q + (5 n ) , (49) 

3=1 

where in the case / > the constant C\ can be found explicitly: 
Q = 2 r(Z + 1)(2tt)-'- 1 C(/ + 1) cos I^i±!2 

(here £(•) is the Riemann zeta function). The preceding estimates imply that 

n = K5; i r - 1 T(r- + l) + 0(l). (50) 
from which it follows that 

*• = {kWTT } +0W 

l/(r+l) 



24 



and that 



e ncr " ~ 



y n exp (n r/{r+1 \KT{r + l)) 1 /^)) . 
The asymptotic for B\ follows from 

J' = l ^ ' 

^ (\ — y-j e -js n y _ y-j e -js n y 

n 

~ ^r(r + 2)5 T ; r - 2 

^ K-V(-+l) r ( r + ^"(r+^+llp^ + 2)n ( r+2 )/(r+l) 

The second factor in (14) may be expanded as 

n / n 

] J (1 - e' 1<T ") ~ aj = exp I - a i Ml - y~ j e~ jSn ) 

3=1 \j=l 



= exp 



a 3 y jk e jk5n 



EE 



k 

d=i k=i 



= exp J] «iir'e-"* + EE ) (51) 

\j=l j=l k=2 J 

We use (49) and (50) to show that the first term in the exponential in (51) 
equals 

n n n 

j=l j=l 3=1 



oo 
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where 



D r = iT 1 /( r+1 )r(r)(r(r + l))"'/^ 1 ) = -(KT(r + l) 1/(r+1) . 



The second term in the exponential in (51) equals 



3=1 k=2 ' 3=1 k=2 



where the double sum on the right converges absolutely because of (48). 



5 Asymptotics for expansive selections 

Let c* be the number of selections of total size n corresponding to a given 
sequence dj. The generating function for c* is given by 



?1=0 

OO 



n=l ffen n j=l x '■> 



3 \x j v 



= i[(i+x j ) aj ,\x\ < 1. 

By adapting the derivation of (4) for multisets to the truncated generating 
function g* n (x) = YYj=i (1 + we obtain for all a G R, 

c* n = e na / ]"T (1 + e - ja e 27viaj ) aj e~ 2man da. (52) 
Jo j=1 

It follows that 



n pi n / 1 1 p -ja p 2iriaj \ <*j 

n (1 + I n ( ) e " 2 """< i « 

71 „1 

= e ncT TT (1 + e~ ]a ) aj / 0*(a)e- 2 ™ n rfa, 
.7=1 Jo 



(53) 
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where 

n 

(f>*(a) = Y[(f)*(a), ae 

3=1 

and 

/ 1 _|_ p -j<r p 2maj \ "3 



1 + 

If a > 0, then the 0* are characteristic functions of a sequence of independent 
binomial random variables j -1 ^- : 



r(X k = JJ) = (7) ' ( TT L C ) ■"',/ = 0, 1, 2, ... , a, 

The formula (53) could also be derived from (145) of [3]. 

The number of integer partitions of n with distinct parts all of size at 
least s was considered in [13]. This is the selection with 



a j 



if j < s, 

1 if j > s. 



The identity (53) was derived in [13] for this particular example. 
Let Y = YT i=1 Xj. We have 



m:-.= ey = j2{ 



n . _ ■ 



We will assume that a = a n is chosen in such a way that M* = n. The 
fact that a n can be chosen in such a way follows from considering that M* 
decreases from \ YTj=i 3 a j to 0> as ° changes from to +oo and noting that 
the assumption that the dj satisfy (13) implies that YTj=i3 a 3 > n for n large 
enough. Under the above choice of a, the variance of Y is 

? *\2 V / V \ _ J U 3 C 



3=1 



From this starting point the proof of Theorem 5 is similar to the proofs 
of Theorem 1 and Theorem 4. 
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